Signal processing using wavelet transform for MEG noise reduction 
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1 Introduction 

The noise involved in data gives large uncertainty 
on the source localization of MEG. The signal to 
noise ratio (S/N) of MEG data evoked by olfactory 
or taste stimulations can not be improved by 
averaging because of subject’s habituation. To 
increase the accuracy of the source localization, we 
developed a noise reduction method using wavelet 
transform, because the wavelet transform was 
suitable for time-frequency analysis. 

We assessed this noise reduction method using the 
real auditory evoked MEG (AEF). The reason is as 
follows. The source location of the AEF is 
anatomically obvious and we can obtain the data 
with any S/N by changing the number of averages. 
The source locations with only 5-time averaged AEF 
data (S/N was nearly 2) were accurately estimated 
by applying the noise reduction method and the 
estimated positions were agreed with the anatomical 
position. 

2 Methods 

2.1 Wavelet transform and mother function 

FFT (Fast Fourier Transform) is widely used for 
signal processing. However, it is not effective to 
noise reduction of MEG data. Because, the sine and 
cosine, which comprise the basis of Fourier analysis, 
are non-local in time space and the temporal 
information is lost by the transform. 

Toward this, the functions that have compact 
support (mother wavelet functions) are used as a 
basis function of wavelet transform. An advantage 
of the wavelet transform is that the basis functions 
can be varied to suit one’s object. 

We adopted the “Daubechies 8” is the basis function 
as a mother wavelet function. The scaling function is 
the basis function for the mother wavelet. The 
scaling function plays as the low-pass filter and 
mother wavelet function plays as the high-pass filter. 
Figure 1 shows the “Daubechies 8” mother wavelet 
function. The frequency-amplitude character of 
Daubechies mother wavelet family is shown by the 


form of the Daubechies N (N=l,2,3...) in Fig.2. As 
the N becomes larger, the characteristic of frequency 
becomes better. 



Figure 1: (( Daubechies 8 ” mother wavelet function 
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Figure 2: Frequency-amplitude characteristic of 
‘Daubechies 8” 

2.2 Multiple resolution analysis 

We adopted the multiple resolution analysis with 
wavelet transform for reduction of noise [1, 2]. In 
the multiple resolution analysis, the data set was 
decomposed to sub sets with each frequency 
resolution. They have a pyramidal structure as 
shown in Fig.3. The suffix n presents the level of 
frequency scale. 

The frequency component c 0 is decomposed to Ci and 
di. This process is repeated in order. That is, d n is 
decomposed to c n +i and d n +i, where c n +i is the lower 
frequency component of c n and d n +i is the higher 
frequency component of c n . The original signal is 
represented in terms of all c n and d n . Decomposition 
and re-building are given with the following 
equations. 





Figure 3: Multiple resolution analysis 
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standard deviation of the frequency distribution was 
used as the threshold A, 5 [3]. 

a) Auditory evoked MEG data 
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b) High frequency element d 5 of the data 
shown in Fig. 4(a) 
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where, 

c 0 : original data 

c n : the low frequency component of the level n 
d n : the high frequency component of the level n 
p n : the low-pass filter of level n 
q n : the high-pass filter of level n. 


2.3 Thresholding 
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c) Frequency distribution of d 5 
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Noise reduction is carried out by substituting zero 
for high frequency components (d n ). If we substitute 
zero for all high frequency components, the signal 
components included in d n are neglected. To re¬ 
build the more accurate signal, “thresholding” was 
used. The idea of thresholding is to substitute zero 
for all elements in d n , d n k that are less than a 
particular threshold A". 
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Since the A n should be determined without 
arbitrariness in each process, we decided A” as 
follows. Figure 4 (a) shows a MEG data evoked by 
auditory stimulation. The high frequency component 
d 5 of the MEG data is shown in Fig. 4(b). The 
frequency distribution of the amplitude shown in Fig. 
4(b) is shown in Fig. 4(c). This figure shows that the 
elements in the d 5 separate into two; groups, one 
group is near center and has a distribution similar to 
the Gauss distribution and the other is apart from the 
former group. We can recognize the elements in the 
former group as noise elements and the elements in 
later group as signal elements. Considering this, the 


Figure 4: Decision of threshold A 5 : a) Auditory 
evoked MEG data, b) High frequency element d 5 of 
the data shown in Fig. 4 a), c) Frequency 
distribution of d 5 . 

3 Application to auditory evoked MEG 

3.1 AEF 

To examine the validity of this noise reduction 
method, we measured the auditory evoked MEG at 
Electrotechnical Fab., Osaka, Japan by Neuromag- 
122. The subject was a 25-years-old and healthy 
male. The auditory stimuli of tone burst of the 
frequency 1000 Hz were applied to both subject’s 
ears. 


3.2 Noise reduction 

5 data; not averaged data, 5 times, 10 times, 50 
times and 100 times averaged data, were prepared 
for the analyses from the measured MEG data as 
original data set for source localization. Figure 5 
shows 5 times averaged AEFs on the 122 channels. 
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Figure 5: Auditory evoked MEG after five times 
averages 

The averaged AEFs were decomposed into level 5, 
were applied the thresholding and were re-built. The 
re-built AEFs are shown in Fig. 6. A peak clearly 
appears at the latency of about 100ms. 



Figure 6: AEF after the noise reduction 

The S/N is a kind of indicator on the effectiveness of 
the noise reduction method. The changes of S/N due 
to the noise reduction by this method are shown in 
Table 1 as inverse number of S/N, where we define 
N/S by the following equation. 

(The standarddeviationof the noise)x 3 

N/S = --—x 100[%] 

| The maximumamplitudeof the signal | 


Table 1: Improvement of N/S 


number 

of 

averaging 

N/S before noise 
reduction [%] 

N/S after noise 
reduction [%] 

0 

104.38 

74.09 

5 

47.06 

31.70 

10 

36.45 

26.80 

50 

18.89 

16.09 

100 

18.04 

17.02 


3.3 Source localization 

We estimated signal sources to evaluate the validity 
of the presented method for noise reduction. The 5 
data sets; not averaged, 5 times, 10 times, 50 times 
and 100 times averaged data, were prepared as 
original data. For comparison, the source 
localization was done on these 5 data sets before and 
after noise reduction by the presented method. 


A combined method of GA and SA was used for the 
source localization with 2-dipole and spherical 
model [3]. The estimation was carried out in the 
period of 20 ms (2.5 ms interval, 9 points) around 
the latency of 100ms. 

3.4 Result 

The locations and the moment directions of the 
estimated dipoles with 5 times averaged original 
data and with the data after noise reduction were 
shown in Fig. 7 and Fig. 8 respectively. The 
positions scatter within 35 mm in the case of 
original data, though the scattering range is within 
several mm in the case of noise reduction data. 
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Figure 7: Result of the 
data. 
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Figure 8: Result oj the 5 times averaged data with 
noise reduction. 


4 Discussion 

Table 2 shows the scattering range of the estimated 
position with each data. The scattering range 
becomes small by applying this method. Especially, 
the efficiency of the method is high on low number 
averaged data. Because, the high number averaged 
original data itself did not include so much noise. 
The estimated position after noise reduction is 
agreed with the center of the scattering range for the 
original data as indicated in Fig.7 and Fig.8. The 
estimated dipoles after noise reduction are overlayed 
on his MRI and shown in Fig.9. The estimated 
dipoles are nearly on the auditory area. These results 
show this method is useful for source localization on 
the MEG with low signal to noise ratio. 





Many times averaging required much measuring 
time and there is a risk that the subject’s posture 
changes during the time. The posture change causes 
some shits in the estimated position on the many 
times averaged MEG data. The presented method 
may be effective to avoid the error in source 
localization. 


Table 2. Scattering range of the estimated position 


Average 

number 

Before the noise 
reduction [mm] 

After the noise 
reduction [mm] 

Right 

hemisphere 

Left 

hemisphere 

Right 

hemisphere 

Left 

hemisphere 

0 

43.57 

23.50 

23.52 

11.64 

5 

25.81 

21.18 

7.44 

6.56 

10 

7.81 

20.83 

5.95 

13.04 

50 

6.67 

9.52 

4.81 

9.57 

100 

4.44 

8.35 

3.76 

9.28 


a) dipoles in the left himesphere 



5 Summary 

To increase the accuracy of the source localization 
on MEG with low signal-to-noise ratio, such as 
olfactory or taste MEG, we developed a noise 
reduction method using wavelet transform. 

This method was assessed using the real auditory 
evoked MEG. We prepared 5 original data that had 
different S/N by changing the number of averages. 
The noise reduction rate and the change of estimated 
error with source location were investigated by 
application of this method. It is made clear that the 
accuracy of source localization is improved by 
applying this method. 
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Figure 9: Overlay of estimated dipoles on his MRI 
a) dipoles in the left himesphere b) dipoles in the 
right hemisphere. 

















